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The collision smearing of the nucleon momenta about their initial values during relativistic nucleus-nucleus 
collisions is investigated. To a certain degree, our model belongs to the transport type, and we investigate 
the evolution of the nucleon system created at a nucleus-nucleus collision. However, we parameterize this 
development by the number of collisions of every particle during evolution rather than by the time variable. 
It is assumed that the group of nucleons which leave the system after the same number of collisions can be 
joined in a particular statistical ensemble. The nucleon nonequilibrium distribution functions, which depend on 
a certain number of collisions of a nucleon before freeze-out, are derived. 
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1. Introduction 

The problem of isotropization and thermalization in the course of coUisions between heavy relativis- 
tic ions attracts much attention, because the application of thermodynamic models is one of the basic phe- 
nomenological approaches to the description of experimental data. Moreover, the assumption regarding 
a local thermodynamic equilibrium, along with other factors, is successfully used in various domains of 
high-energy physics. Meanwhile, many questions concerning this problem remain open for discussion. 

The main goal in the investigations of the collisions of relativistic nuclei is to extract the pertinent 
physical information on the nuclear matter and its constituents. In the present paper we propose a trans- 
parent analytical model of the nucleon spectrum which occurs in the course of relativistic heavy-ion 
collisions. Our model is aimed at extracting the physical information from the nucleon spectra which 
concerns the nonequilibrium processes. 

Let us look at the cross-section of a multiparticle production during the collision of two nuclei "A" 
and "B" (see figurejl). In order to describe of the nucleon subsystem one can parameterize the final state 
of the nucleon ensemble by 4-vector - (^w, ^^jv)- In the center of mass of this Af-nucleon ensemble, 
where - 0, the total cross-section reads 



div 
O'nucleon 



JV 



d3pid3p2 ■■■d^pjv „=i 



n ^(Pn) 



= W(p^,P2,...,Pn)S 



En 



N 



n=l 



(1) 



where = x/m^Tp^ (in the final state, the particles are regarded as free ones) and we adopt the 
system of units where the speed of light c - I. Due to the presence of the 5-function, which "fixes" the 
energy of the nucleon system, the last expression (ij looks like a probability in the microcanonical ensem- 
ble. Then, it is reasonable to make the Laplace transform with respect to the energy of the nucleon 
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It turns out that now one can describe the final state of the nucleon subsystem through one of the two 
parameters: the total energy Bjv or the parameter p. 

All the above formulae were introduced for a brief discussion of the basics of the statistical model (for 
details see fT/ZI). Actually, the statistical model neglects all the dynamics accumulated in the transition 
probability W in favor of the features of the phase space. Formally this is expressed like approximation: 



W[p^, P2,---, Pn) ~ const. 
Then, from j2} one immediately obtains the multi-nucleon cross-section parameterized by /3: 

d O^nucleon ^ ^ 



d^pid^P2 ...d^PN „=i 



N 

n ^(Pn' 



n 



(3) 



(4) 



where ¥/„ are some constants and Yl^^i Wn - const. 

On the other hand, the statistical description, which arises after freeze-out, is conceptually restricted 
just to several conserved quantities: total energy of the system En, number of particles N, and conserved 
charges such as the baryon number. Of course, this dependence can be transferred to the descriptions 
by means of conjugate quantities: parameter p and chemical potentials which are in correspondence 
with N and the conserved charges. As we see from j4f , this description provides a certain information on 
the spectrum of the registrated particles. Meanwhile, any dynamical information on multiscattering pro- 
cesses during collisions is lost. At the same time, it is well understood that a microscopic description can 
be carried out just on some level of approximation. For instance, if it is possible to factorize the transition 
probability W, i.e., to write it in the form Wip^, p2, p^) ~ X\^=\ W[p^) , we come to factorization of 
the multi-nucleon cross-section for the particles of one species 



d^Pld^P2-..d^PN n=i 



N , , 

= n [mp„)e-^'""'"']. 



(5) 



n=l 



The approximation of the sequential rescatterings of a particle during collision of nuclei which is pro- 
posed in the present paper is exactly in this framework. We follow the chain of reactions (rescatterings) 
of every nucleon that goes through a number of hadron transformations, and we obtain a single-particle 
spectrum of the nucleon (p) which depends on the number M of collisions (reactions) of the nucleon 
(of the baryon characterized by the baryon number B - 1). 

We argue that the number of nucleon collisions (reactions) at AGS and SPS energies is finite and the 
maximal number of collisions Mmax is not more than Mmax = 43. Apart from this, all the nucleons which 
are freezed out during a particular nucleus-nucleus collision, can be subdivided into groups. In every 
group, the nucleons went through the same number M of collisions. We determined that the most popula- 
ted groups are for the number of collisions M which are in the range: M = 14 - 17. Starting from the initial 
moment of the nucleus-nucleus collision, we follow the sequential collisions of every nucleon through the 
opposite nucleus (see figure[lj. Nevertheless, the original nucleon can be transformed during a particular 
collision into another particle, for instance into delta isobar . Then, we follow a new particle which 
carries the same charges (the baryon number, electrical charge, etc.) as the original nucleon. During the 
last collision (it can be a decay), all these "intermediate" particles transform back into nucleons. Hence, 
we can investigate just nucleons in the final state. Starting from this point, every group of nucleons in the 
final noninteracting state is the subject of a statistical model. For the multiscattering stage of evolution, 
we treat UrQMD l3j|4J as a source of "experimental data" which we use to adjust the parameters of our 
model. 

We propose a mutually complementary combination of these two approaches, i.e., an approximate 
description of the dynamical stage of evolution of nucleons during the nucleus-nucleus collision which is 
completed with a statistical description of the nucleon freezed out stage. Our approach is based on "The 
multiscattering-statistical model" elaborated by us. 
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2. The multiscattering-statistical model 

Consider successive variations of the momentum of a nucleon from nucleus A (see figure [ij which 
moves along the collision axis from left to right through the nucleus B. Every m-th collision induces the 
momentum transfer, q^, for this nucleon. Consequently, after M collisions, the nucleon acquires the 
momentum k: 

fco — fco + ^i ko + q-^ + q2 — ■■■ — ko + Q-k. (6) 

where Q - E^^i 1m ^^e total momentum transfer finally obtained by our nucleon after M collisions 
(see figure [2|. If the M-th collision is the last one, then the nucleon is free after having been freezed out 
from the system. As a matter of fact, it will be a group of such nucleons which experienced the same num- 
ber of collisions M before the freeze-out. The relation of these nucleon groups (nucleon sub-ensembles) 
to the spectrum is discussed in the next section. 




Figure 1. Two colliding identical nuclei. Two- 
stream system is created during the collision of 
every nucleon from nucleus "A" with nucleons 
from nucleus "B" and vice versa. 



Figure 2. Transformation of the initial nucleon mo- 
mentum, fco, as a result of M collisions; Q = 
^m=i '/m is '^o'^sl momentum transfer after M 
collisions, is the momentum transfer in the m-th 
collision. 



We assume that all the momentum transfer obtained by the nucleon from the nucleus "A" during 
its travel through the system "B" are some random quantities. Then, when the number of collisions M 
is big enough in accordance with the central limit theorem, the distribution of the random quantity Q 
should be a Gaussian one 



GiQ)-- exp 
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{Q-(Q)f 

20-2 



1 



- exp 



(fc-fco-<Q>Mf 



2a 
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(7) 



where we take into account the equation |6|, i.e., Q = fc - fco, and explicitly write the dependence on 
the final momentum of the nucleon fc. Here, z is the normaUzation constant. Actually, this distribution 
depends on the number M of random quantities q^, q2, 1m' which coincides with the number of 
collisions M experienced by the nucleon before being freezed out. 

Further, we design the many-particle distribution function Gm(£'m; fci , fc2, ■ ■ ■ , fcjv) as 



Gm{Em; fc) = 



1 



k- (fci, fc 



2, ... 
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and make the Laplace transform with respect to the total energy Em 



(8) 



(9) 



(10) 
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Basically, from now on, any physical quantity that depends on the nucleon momenta, can be averaged 
using the many-particle distribution function 

GmW; h = — ^ GMiP; h . (11) 

If we introduce expression j9j into (lOj, then in the framework of the multiscattering-statistical model 
(MSS-model) we obtain 

1 rT„„„l « .ui. . [kn-ko-(Q)M) 
where Z^ifi^) - [z„iPJ^,)]^" and the single-particle distribution function is as follows: 



Gm(j6m; = ^TTT^ n exp \ -P^c^iK) - ^ ^2 \ = U f^^^n). (12) 

^m^Pm> n=l \ '^'^M 



1 fc-fco-<Q>M) 

/M(fc) = — ^ exp \ -P^co(k) - ^ " , \ (13) 

with Zj^^ ()3j^) as the single-particle partition function. Note, a derivation of the analogous distribution can 
be found in (Dill. 



3. Two-stream dynamics 

Based on the obtained results, we are coming to a description of a two-stream system which is cre- 
ated in the course of relativistic nucleus-nucleus collisions. The description is based on the following 
assumptions: 

1. Starting from the initial state (first touch of the colliding nuclei), at an arbitrary moment of time, 
there are two systems of nucleons: one system consists of nucleons with a positive z-component of 
the nucleon momentum, i.e., Pz S (we refer to this system as "A") and the second system consists 
of the nucleons with a negative z-component of the nucleon momentum, i.e. Pz^O (we refer to this 
system as "B", see figure[lj. Even after the freeze-out, there is a good enough separation of these 
systems along p^-axis. 

2. The number of collisions of every nucleon (hadron) is finite because the lifetime of the fireball is 
limited. To determine the maximal number of collisions, Mmax, in a particular experiment we use 
the results of UrQMD simulations (3]|D. 

3. Since the colliding nuclei are spatially restricted, different nucleons experience a different number 
of collisions, and it is intuitively clear that the collision histories of the inner nucleons and sur- 
face nucleons will be different. That is why, we subdivide all the freezed out nucleons (nucleon 
ensemble) into different nucleon sub-ensembles in accordance with the number of collisions M 
the nucleons went through before being freezed out. Then, the nucleons from every sub-ensemble 
give their own contribution to the total nucleon spectrum. 

It is time to write a two-stream distribution function Fjvf(p) which, in accordance with the first as- 
sumption, should take into account the spectrum produced from both particle streams, "A" and "B". Being 
normalized to unity, a two-stream distribution function looks as follows: 

^Mip) ^ l[f^\p) + f^hp)], (14) 



where 



and 



Mr,. 1 f «... I 1 {kz-koz-(Qz)M)\ 



.W(fc)-^_exp|-^ !cxp! (^- + ^°-+<Q->Mf 
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with ko = (0, 0, koz) and we assume ((t^)m - (cr^)M - ((^\)m ■ Here, the single-particle partition function 
reads Zj^ip^) - f d^k/(2n)'^f^\k). Now, the slope parameter reflects also a collective motion of the 
M-th nucleon sub-ensemble moving in the laboratory system. 



Au+Au, 10.8, 40, 80 AGeV, Pb+Pb, 158 AGeV 



Au+Au, 2, 10.8, 20, 40, 80 AGeV; Pb+Pb, 158 AGeV 
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Figure 3. (Color online) The populations Nm of the 
nucleon sub-ensembles which depend on the num- 
ber of collisions M. The result is obtained from the 
UrQMD simulations for the most central collisions. 



Figure 4. (Color online) The z-component of the 
mean momentum transfer versus the number of 
coUisions. 



Keeping in hands the two-stream distribution functions Fjvf(p), where M - 1, 2, ... Mmax. one can 
construct the nucleon spectrum which occurs in the course of a central nucleus-nucleus collision. If we 
denote the number of particles in a particular sub-ensemble, where the nucleons experienced M colli- 
sions, by Nm, then in the c.m.s. of the colliding nuclei, the total nucleon spectrum is as follows: 

= E ^mFm(p), with ^ ^M=Ntotai, (17) 

^ P M=l M=l 

where ATtotai is the total number of net nucleons. The number of nucleons Nm in every sub-ensemble cal- 
culated for different energies with the help of the microscopic transport model UrQMD 1 3 4| is depicted 
in figure [3] Every sub-ensemble of nucleons can be described as an ideal gas which moves with some 
collective velocity. 



4. Extraction of physical parameters from the data 

First we obtain from UrQMD f3^ '41 the longitudinal distribution of nucleons for every M-th sub- 
ensemble. The distributions for the stream "A" (positive pz) and the stream "B" (negative pz) were ob- 
tained separately. We refer to this pool of distributions as "UrQMD data". We fit the "UrQMD data" on the 
longitudinal distribution of nucleons of the M-th sub-ensemble of the stream "A" exploiting the theoret- 
ical distribution function (Ts) integrated over a transverse momentum. The variations of the theoretical 
distribution function were provided by four parameters: {Qz)m> Pm> {^\.]m ^^'^ i^^)M' '^^^ results of 
the fit of the "UrQMD data" (nucleon longitudinal distributions) for the energies 10.8 and 158 AGeV are 
depicted in figure[5] 

The dependence of the parameters with respect to the number of reactions, M, of the nucleon before 
the freeze-out is shown in the following figures: The slope parameter T,^ = l//3j^ in figure [6] (the curves 
marked with square symbols) and the temperature Tq (the curves marked with circle symbols); The mean 
value of the shift of the distribution function, {Qz)m> in figure [i] The longitudinal variance in 
figurejs] Note, we use the system of units where the Boltzmann constant is unit, k^-l. 

Temperature of the hot ideal gas Tq is determined in the local rest frame of the gas. This temperature 
is connected with the total kinetic energy of N particles (nucleons) in the local rest frame in the following 
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Au+Au, 10.8 AGeV Pb+Pb, 158 AGeV 




, GeV/c , GeV/c 



Figure 5. The spectrum of the M-th nucleon sub-ensemble (M = 6, 10, 18, 30) with respect to pz- 
component of the nucleon momentum calculated using the UrQMD transport model for Au+Au collision 
(black squares). The grey curves are the fits of the UrQMD data within the framework of the proposed 
multiscattering-statistical model (MMS-model). 



way: 



= 3 To + m„ isTi 

N [To 



To 



(18) 



In our model is the total kinetic energy of the sub-ensemble of nucleons in the rest frame of this 
group of particles, is the nucleon mass (E^j - y[E- V Pz] , j- 1/ Vl- V^). 
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Figure 6. (Color online) Dependence of the slope parameter (curve marked with squares) and tem- 
perature To (curve marked with circles) on the collision number M for different energies of the nucleus- 
nucleus coUision: 10.8, 158 AGeV. 



To obtain the transverse distribution we integrate the "A"-stream distribution function (Ts) over the 
longitudinal component of the nucleon momentum. The results of the description of the "UrQMD data" 
on the nucleon transverse distribution for the energy 20 AGeV is depicted in figure [7] (left hand panel). 
For this description, we use the same values of parameters {Qz)m> Pm> (^z)m ^^'^ i^\)M' which were 
obtained during the fit to the "UrQMD data" on nucleon longitudinal distributions. The experimental data 
for transverse nucleon distributions in central Au-i-Au reactions at pi^t, - 11. 6A GeV/c |7| were described 
using the formula (Tt) , the result is depicted in figure [7] (right hand panel). We see a good agreement of 
the description with experiment. 
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5. Discussion and conclusions 

The description of a many-particle system, which is in thermal equilibrium state, can be regarded 
as quite understandable and complete by means of the temperature and chemical potential if the latter 
is needed. Then, to obtain the value of the temperature, one has to fit the particle spectrum data using 
one of single-particle distribution functions. The fitting procedure is nothing more as an extraction of the 
physical quantity, i.e., "temperature", from the data. 



Au+Au, E=10.8 AGeV Pb+Pb, E=158 AGeV 




collisions M for collision energies 10.8 and 158 AGeV. 



Next, if we come to the description of a many-particle system in a nonequilibrium (nonstationary, 
nonhomogeneous) state, a natural question arises: which set of parameters is needed to get a relevant 
physical picture of the many-particle system, which state evolves in time? Of course, by the words "phys- 
ical picture" one means a physical interpretation of the parameters of a model. The set of parameters of 
a model, as well as the behavior of the evolution of the parameters take the form of a specific language 
based on which we investigate, describe and "understand" our nonequilibrium many-particle system and 
the processes inside it. 

In the present paper for description of a nonequilibrium state we propose three parameters, which 
are defined in "The Multiscattering-Statistical Model". 

1. The slope parameter T(M] - I/^m, which reflects as well a collective motion of the nucleon sub- 
ensemble. Its dependence on the number of collisions M experienced by nucleons, which belong to 



13201-7 



D. Anchishkin, V. Naboka, J. Cleymans 



M-th sub-ensemble, is depicted in figures [6] By means of the energy per particle (after the freeze- 
out), this parameter is directly connected with the instant temperature To{M) of the M-th nucleon 
sub-ensemble which is defined in the frame where the sub-ensemble is in rest. 

2. The mean value of z-component of the total momentum transfer (Qz)m - Zm=i^'?z^'«> which is 
related to the kinematics of a particle-particle collision. 

3. The variance of the Gaussian distribution. Actually due to the different collision dynamics along 
the different axis, the variance is split in two pieces: a) The longitudinal variance [a^]^^. b) The 
transverse variance (c^)^. In figuresjij the behavior of the transverse variance is confronted with 
the longitudinal variance. We see a steady broadening of the transverse distribution with an in- 
crease of the number of collisions M in the nucleon sub-ensemble. Indeed, the initial transverse 
variance of nucleons is approximately zero on the scale of collision energy. The distribution of nu- 
cleons around the initial value in the transverse direction becomes broader and broader with time. 

We can conclude that the evolution of the physical parameters of "The Multiscattering-Statistical 
Model" elaborated in the present paper gives a transparent insight into the dynamics of the net nucle- 
ons in the course of relativistic nucleus-nucleus collisions. The appearance of the Gaussian distribution 
as a factor in the nonequilibrium distribution function of nucleons, see (Ts) , is common in describing any 
multiscattering process with a big but finite number of rescatterings of the particle when we can regard 
every particular scattering independent of others. This condition is especially satisfied for the nucleons 
in the course of high energy heavy-ion collisions when the particle wavelength A - hi pis much smaller 
than the mean distance between the nucleons in a nucleus. 
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HepiBHOBa>KHi ())yHKL4M posnoAi/iy HyicnoHiB npn 

pe/UtTMBiCTimHlIX 5IAP0-5IAepHI/IX SITKHeHHJIX 

P,. AHHi/iLUKii-R B. Ha6oKcP, X. K;iewMaHcP 

iHCTWTyT TeopeTWHHoT c()i3HKW Im. M.M. BoroyiK)6oBa, 03680 Khib, VKpaiHa 
Kwi'BCbKwi^ HaL4iOHa/ibHi/iti yHisepcwTeT iMeni Tapaca LLIeBi-ieHKa, 03022 Khib, VKpaiHa 
VHiBepcMTeT Kewn Tayna, PoHAe6oLU 7701, HiBfleHHO-AcjjpwKaHCbKa Pecny6/iiKa 

floc/iiflxycTbCfl posMHTTfl iMnyyibciB HyK/iOHiB HaBKO/io CBOi'x nonaiKOBnx SHaneHb, ^Ke Bifl6yBaeTbCfl b pe/in- 

THBiCTH4HI/lX flflpO-flflepHMX SilKHeHHaX. HaiUa MOfle/lb BiflHOCHTbCfl, neBHOK) MipOK), flO TpaHCnOpTHMX, MH p,o- 

cn\/i,v\m eBO;iK)i;iK) HyKJiOHHoT cwcreMW, CTBopeHoT b nflpo-flflepHnx siTKHeHHflx, a;ie mi/i napaMerpi/isycMO i4eM 
po3Bi/iTOK He nacoM, a 4i/ic;iom aiiKHeHb koxhoT 4acTHHKW. npunycKaeibCfl, iu,o rpyna HyK/iOHiB, hkI 3a;inujaK3Tb 
CMCreMy sasHaBiui/i o^HaKOBy Ki/ibKicTb siTKHeHb, MOxyTb 6yTH o6'eAHaHi b neBHuti CTarwcTUHHuti aHcaM6;ib. 
06paxoByeTbCfl HepiBHOBa>KHa (|)yHKL4ifl posnofli/iy HyK/iOHiB b iMny/ibCHOMy npocTopi, ^Ka sa/iexnTb Bifl neB- 
Horo 4MC;ia siiKHeHb HyK/iOHa nepefl BwnpoMiHKJBaHHaM is ci/icreMH. 

KniOMOBi cjiOBa: pejian/iBiCTH'-iHi siTKHeHHfi, HepiBHOBa>KHa ipyHKU,ifi posnofliny, cneKTp HyKJiOHiB, 
napaMBTpusaLiifi eBOJiKJLfi'f 
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